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Abstract. Results from a cross-correlation analysis be- 
tween the COBE DMR 4 year, and ROSAT PSPC All- 
Sky Survey data are presented. Statistical comparisons 
between microwave and X-ray maps can probe interest- 
ing astrophysical environments and processes, such as the 
warm interstellar medium, the Sunyaev-Zel'dovich effect 
in clusters of galaxies or gaseous group halos, X-ray lu- 
minous radio sources and the Integrated Sachs- Wolfe or 
the Rees-Sciama effect. In order to test the diffuse, extra- 
galactic X-ray background as probed by ROSAT, against 
the COBE DMR large-scale CMB structure, our analysis 
was performed in most detail in a ROSAT selected re- 
gion of the sky (+40° < b, 70° < 1 < 250°) in an X-ray 
energy range with minimal Galactic structure and resid- 
ual X-ray contamination, and the COBE low noise and 
least Galactic contribution channels. Comparing to other 
regions of the sky and neighbouring energies and frequen- 
cies, we find indication for a positive Galactic correla- 
tion on large scales. This correlation is most prominent 
in the softest X-ray band and lowest microwave chan- 
nel, with a > 95 % confidence level detection against 
COBE noise and CMB cosmic variance, including the 
high quadrupole value resulting from the power spectrum 
fit. The spectral dependences are consistent with Galactic 
thermal X-ray emission, and Galactic synchrotron radia- 
tion or free-free (Bremsstrahlung) emission by the warm 
interstellar medium in the microwave regime. Removing 
the quadrupole term on a sky map with a Galactic cutout 
or related gradients in the selected regions leaves no cor- 
relations above a X-a level on smaller angular scales. We 
conclude that there is no significant extragalactic correla- 
tion on scales for which the combined data are sensitive 
(7° — 40°) and that Galactic correlation is significant only 
on large angular scales, of the order of the quadrupole. In 
the context of removing large angular scale gradients we 



Send offprint requests to: R. Kneissl, 
ruk@mpa-garching.mpg.de 



give results on best fit X-ray dipoles from various ROSAT 
data and discuss these with respect to the difficulty of 
finding a cosmological dipole. The lowest correlation up- 
per limits we can place are ~ 15% of COBE CMB fluc- 
tuations and ~ 10% of the ROSAT extragalactic XRB 
fluctuations. We discuss these results with respect to the 
possible correlation mechanisms. 

Key words: galaxies: clusters: general — cosmology: cos- 
mic microwave background - cosmology: observations - 
cosmology: diffuse radiation — X-rays: general 



1. Introduction 

Astrophysical radiation backgrounds are an important 
probe in cosmology. The instrument data used in this anal- 
ysis, taken by the Cosmic Background Explorer's (COBE) 
Differential Microwave Radiometer (DMR) (Smoot et al. 
1990) and the Rontgen Satellite's (ROSAT) (Triimper et 
al. 1990) Position Sensitive Proportional Counter (PSPC) 
(Pfeffermann et al. 1986) All-Sky Survey (RASS) have 
contributed to the improvement of our knowledge of the 
cosmic microwave (CMB) and the cosmic X-ray back- 
ground (XRB). The COBE DMR discovered for the first 
time structure of cosmic origin in the CMB (Smoot et al. 
1992). The ROSAT PSPC data helped to resolve a large 
fraction (more than 75%) of the emissivity of the XRB 
into point sources with redshifts out to z ~ 3 (Hasinger 
et al. 1993, Comastri et al. 1995). 

As the two experiments probe very different energy 
regimes (Tcmb = 2.7 K, Tik e v — 10 7 K, Planck tempera- 
tures) and the backgrounds originate at different redshifts 
{zcmb ~ 1000, zxrb ~ 0-3) no strong correlations can 
be expected between the two data sets. However physical 
processes connecting the microwave and X-ray radiation 
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and common emitting sources do exist, leaving the possi- 
bility for a weak correlation. 

Both data sets show the Galactic plane prominently 
and even at high Galactic latitude Galactic emission is 
present in the maps (Kogut et al. 1996). Correlation in this 
case can mean spatial alignment of the emitting regions 
or common physical processes. Galactic contributions to 
the microwave radiation in the DMR frequencies are syn- 
chrotron radiation in the Galactic magnetic field, dust ra- 
diation at 18 K, and Brcmsstrahlung from the warm inter- 
stellar medium of 8000 K (Bennett et al. 1992). The diffuse 
Galactic X-ray radiation consists of continuum and line 
emission from hot plasma of 1 to several million K. Struc- 
tures are mainly constituted by the emission of nearby su- 
pernova remnants and superbubblcs in the harder bands 
(~ 1 kcV) but also from absorption through neutral gas 
associated with Hi clouds in the softer bands ( < 0.5 keV) 
(Snowden et al. 1995; Egger et al. 1996). 

An astrophysical environment directly connecting mi- 
crowave and X-ray radiation is a hot plasma with high 
column density along the CMB photon path, as was first 
noticed by Sunyaev & Zel'dovich (1971) (SZ). In clusters 
of galaxies, electrons with temperatures of <~ 10 8 K and 
densities of ~ 3 x 10~ 3 cm~ 3 scatter photons of the CMB 
in the inverse Compton process to higher energies leading 
to y-distortions of the spectrum and a decrease in photon 
number in the Rayleigh-Jeans regime. The same plasma, 
on the other hand, strongly emits at Bremsstrahlung en- 
ergies of 10 keV. Although the dependences on electron 
temperature T e and electron density n e arc different for 
the temperature decrement, AT$z, and the X-ray flux, 
Sx, according to the approximate relations 

AT S z oc n e T e 
_i 

Sx oc n 2 e T e 2 

a strong anticorrelation is present and has been observed 
in individual clusters (e.g. Birkinshaw et al. 1991). 

Recently Suto et al. (1996) suggested a similar effect 
for a possible plasma halo of the Local Group. As the 
assumed temperature (~ 1 keV) and density (<~ 10~ 4 
cm~ 3 ) are lower, the expected effect is smaller than for 
galaxy clusters, but has a predicted angular structure. In 
the COBE DMR data this structure could not be found 
(Banday & Gorski 1996). A combination with the ROSAT 
X-ray template, sensitive in this energy range, further im- 
proves the limits, also testing for a non-spherical structure 
of the halo. 

Since the XRB does consist to a large extent of point 
sources, of which the largest fraction are active galaxies 
(Hasinger et al. 1993), a considerable number are also ra- 
dio sources, e.g. radio-loud AGN. Laurent-Muehleisen et 
al. (1996), for example, find 2,127 sources common to the 
RASS and the Green Bank 5 GHz radio catalog. Correla- 
tion analysis can thus test the fraction and angular scale 
of correlation at which these populations contribute to the 
data. 



Also, recently Turok & Crittenden (1996) proposed 
the Integrated Sachs- Wolfe (1967) (ISW) effect in a non- 
flat cosmological model as a possible source for a cor- 
relation between CMB and XRB data. This effect, for 
the non-linear evolution of a single gravitational poten- 
tial well, also called Rees-Sciama (1968) effect, introduces 
anisotropies into the CMB via the change of the poten- 
tial in time 4> along the photon path in the direction n 
between recombination r rec and reception To time. 

&n) = 2 f° 6(T,n(r -r))rfr. (1) 

Jr rcc 

They calculated the linear evolutionary effect in detail for 
different A-models, with calculations for an open model 
added by Kamionkowski (1996), and found the largest 
contribution at z < 2. Thus, they suggest the distribution 
of sources in the unresolved XRB to be a tracer of the po- 
tential for the redshift interval of interest. If the redshift 
distribution and the biasing of the sources is known, then 
a detected positive correlation can quantify A. A summary 
of all these effects (table 2) and a discussion with respect 
to the results of the correlation analysis is given in section 
5. 

A statistical comparison between observations of the 
CMB and the XRB was first carried out by Boughn & Ja- 
hoda (1993) comparing the 19.2 GHz survey with HEAO-1 
A2 (~10 keV, 3° resolution), and they found no signif- 
icant correlation, based on Monte Carlo simulations for 
noise properties. Bennett et al. (1993) in cross-correlating 
the 1 year DMR data to HEAO-1 found no significant 
correlation for |&| > 30° and with the LMC masked. In 
the 4 year DMR data analysis by Banday et al. (1996) 
an expansion in orthogonal functions on a cut sky and 
a likelihood analysis for the coupling constant between 
the DMR and HEAO-1 data was used in a simultane- 
ous fit to the CMB power spectrum. Again, no signifi- 
cant correlation was found, when applying a specially de- 
signed Galactic cut based on correlations obtained from 
the DIRBE 140 /xm map and masking of the LMC. Using 
the ROSAT PSPC extends previous work to softer X-ray 
energies, higher angular resolution and better sensitivity. 

Additional interest in this analysis arises from the de- 
tection of a spatially extended X-ray source around clus- 
ters of galaxies found in a correlation analysis by Soltan et 
al. (1996a) between Abcll clusters and the ROSAT diffuse 
XRB. Correlating to the COBE DMR can constrain a gas 
halo model for the extended component. 

In section 2 we introduce the COBE and ROSAT data 
sets respectively used for this analysis and explain how 
they were prepared. In section 3 the correlation method is 
described, including our error estimation. Section 4 con- 
tains the results on various angular scales and their de- 
pendences with energy or frequency, and in section 5 we 
discuss the results in the context of theoretical expecta- 
tions of possible correlation mechanisms. 



R. Kneissl et al.: Search for correlations between COBE DMR and ROSAT PSPC All-Sky survey data 



3 



2. The Data 

2.1. COBE DMR 

The COBE DMR measures the sky differentially in 3 
frequency channels (31.5, 53 & 90 GHz). The maps are 
binned in 2.6° x 2.6° pixels, which are considerably smaller 
than the beam width of 7° (FWHM). The CMB signal has 
an amplitude of 35 ± 2/iK and the COBE DMR mean sen- 
sitivity is 26 /iK per resolution element, corresponding to a 
signal to noise ratio of ~ 0.5 on pixel level. The noise level 
varies from channel to channel (31A: 248, 31B: 316, 53A: 
86, 53B: 101, 90A: 146, 90B: 116; in per pixel). The re- 
sults of the COBE 4 year data analysis by the COBE team 
are summarized and referenced in Bennett et al. (1996). 
We used various maps of the COBE DMR data set. The 
final analysis uses the 4 year data, but consistency tests 
with the 1 year and 2 year data were performed. From 
the channels A and B, which are differentially measuring 
the signal, we constructed inverse noise weighted (A+B)/2 
and frequency combined sum maps to minimize the noise. 



ATi, 



( Wi ATi + Wj ATj) 



(2) 



with — l/a* and 1/Wij = l/(w, +Wj), where w* 
is evaluated on the custom cut sky (Banday et al. 1997), 
considering the noise level and exposure. 

All maps used were converted from antenna to Planck 
temperatures. The standard frequency combination for 
our analysis is the 53+90 GHz map, with low noise and 
little Galactic contribution (Kogut et al. 1996). The indi- 
vidual frequency maps and the linearly combined galaxy 
reduced maps (cmb & smb) were used for comparison. 

2.2. ROSAT PSPC All-Sky Survey 

The ROSAT satellite covers a large energy range (0.1-2 
kcV) in the soft X-rays. The harder part is divided into 
4 bands with maximum responses at the following ener- 
gies (R4: 0.7 keV, R5: 0.8 keV, R6: 1.1 keV, R7: 1.5 kcV). 
All bands, particularly the neighbouring ones, have con- 
siderable overlap with each other of up to 50%, due to the 
limited spectral resolution of the proportional counters. 

The RASS intensity (I) and noise (ct) maps were con- 
structed as 



Ct-B 
Ex 



Ct 



a = 



Ex 



(3) 



(4) 



The abbreviations denote count number of received pho- 
tons (Ct), modeled contamination (B) and exposure (Ex). 
The X-ray contamination in ROSAT consists, with vary- 
ing contribution in the different energy bands, of solar 
scattered X-rays, "short- and long-term enhancements" 
and particle background. The noise maps are calculated 



according to Poisson statistics due to photon number 
limitation. The energy band R6 (0.73 - 1.56 keV) is re- 
garded as the best probe for the diffuse cosmological 
XRB, because the systematic uncertainty induced by fore- 
grounds, such as non-cosmic photons, and contamination 
by charged particles is minimized. We concentrate our 
analysis and results to this band, but investigate system- 
atic effects by comparison with the other hard bands. Par- 
ticularly the R5 band is also low in non-cosmic photons 
and in contamination by charged particles, but which con- 
tains, compared to the R6 band, increased Galactic fore- 
ground and can thus hint at a discrimination between the 
Galactic and the extragalactic signal. In spite of the care- 
ful corrections for exposure and elimination of non-cosmic 
backgrounds (Snowden et al. 1995) the final count rate 
distribution is not completely free from residual contami- 
nation. We therefore tested extensively for correlations in- 
duced by exposure or contamination corrections and found 
no negative effects on the results. Unlike the RASS data 
described in Snowden et al. (1995), the data set used here 
has been constructed on a "photon by photon" basis, i.e. 
the intrinsic resolution is that of the detector (<~ 1'). Thus, 
the sensitivity has been improved by avoiding "crosstalk" 
from bright sources. 

The maps used in our work were binned into 0.7° x 
0.7°, with point sources included, in order to compare to 
the complete integrated flux. In a second step they were 
rebinncd to COBE DMR pixel size with varying exclusion 
thresholds for bright sources. The mean intensity of the 
XRB in the R6 is ~ 1.9 cts s -1 per pixel with a fluctuation 
level of ~ 0.17 cts s _1 per pixel and a signal to noise ratio 
of - 2 at COBE DMR pixelization. 

Even the high energy R6 band is, in large regions of 
the sky, dominated by Galactic emission. A field almost 
free from Galactic structures, which is sufficiently large, 
was chosen (+40° < b, 70° < 1 < 250°), hereafter called 
the selected NGP field. This field is the largest simply 
connected patch of the sky probing primarily the XRB. 
Properties of the XRB in this field have been studied in 
a series of papers (Soltan et al. 1996a; Miyaji et al. 1996; 
Soltan et al. 1996b). 



3. Correlation Method 

3.1. Correlation Function 

To minimize the Galactic contribution, a patch of the sky 
which has the highest sensitivity to the diffuse, extragalac- 
tic XRB was chosen. Due to the size (~ 8% of the sky) 
and the peculiar geometry of the patch, a local statistical 
measure for similarities in structure, the 2-point correla- 
tion function, is preferred over global measures such as 
correlated power spectrum components in e.g. a spherical 
harmonic expansion. The form of the correlation function 
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used in this analysis is the Pearson product moment cor- 
relation coefficient C(a) = 



< X{Tj > a — < Xi > a < Tj > c 



y/< X? > a - <Xi >\ 



<Tf> c 



< T > 2 



in an unweighted scheme. Inverse noise variance weight- 
ing has also been used and was found to give unchanged 
results. The correlation coefficient varies between 1 and -1 
for correlation and anticorrelation respectively with in- 
dicating no correlation. For pixel sizes comparable to the 
resolution limit, the statistical uncertainty of the measure 
is (1 — C 2 (a)) ^JN{ij} — 1. The average is taken over all 
pixel pairs {ij} with separation a in the patch. The sub- 
script a denotes that all the terms were evaluated sepa- 
rately for each angle bin. This ensures a correctly weighted 
normalization even in cases when, due to limited area and 
boundary effects, the zero-lag field properties are not a fair 
ensemble average for all angle bins any longer. Note that 
through our choice of the cross-correlation function, the 
bins arc completely statistically independent, and hence 
any apparent correlation between the bins is due to the 
structures in the maps. 

To determine the uncertainties, which we assume to 
be dominated by the DMR noise (section 2) and the cos- 
mic variance of the CMB structure, different techniques 
were applied. We applied a simple method, which intro- 
duces little prejudice (just assumes rotational invariance 
of the data), to correlate to random samples drawn from 
the maps by rotation. Since the ROSAT maps are known 
to contain a strong Galactic, not rotationally invariant 
contribution on major parts of the sky, in contrast to the 
COBE maps, which have been investigated and found to 
be primarily consisting of CMB structures (Kogut et al. 
1996), and since the COBE data predominantly introduce 
the errors, those were rotated around the NPG, and mir- 
ror image rotated around the SGP in 10° steps to produce 
35 and 36 random samples each. The error estimates in- 
duced by this method agree well with our second method 
using simulations of DMR maps. The CMB structure was 
taken to be a random Gaussian field on the sky with a 
power-law (Q rms -ps = 15.3 ^K, n = 1.2) power spectrum 
(Gorski et al. 1996) convolved with the DMR filter func- 
tion (Kneissl & Smoot 1993). The modes used to construct 
the map cover a range in multipole index I = 2-25. Our 
results were compared against <~ 1000 simulations. The 
DMR noise is given as Gaussian pixel noise distributed 
according to the coverage. 



3.2. Spherical Harmonic Fit 

Although we use a correlation function analysis and not 
a power spectrum correlation analysis, it is important to 
study how the correlation results are influenced by par- 



ticularly the low order multipoles. For this we expand the 
sky maps 



+e 



(5) x i^<P)=^2 b em Y em ($,(p) 



(6) 



1=0 m=-l 



into real valued spherical harmonics 



2£+l [ 



sm \m\tp, 
cos \m\ip, 



(l+\m\)\ 1 I 



m < 
= 
> 



p^ix) = (-i)-(i-^) 



dx r 



7 P e (x), m > (7) 



in a simultaneous fit up to multipole order £ max , and sub- 
tract these multipoles from the maps. A well-known prob- 
lem (e.g. Bunn et al. 1994) arises in this procedure. The 
orthogonality relation for the spherical harmonics does not 
hold in the case of incomplete sky coverage: 



L 



(8) 



where Ww m m' ^ o"u'5mm' in general, if R C S\. In our 
case this leads to the fact that we subtract a function 
from the sky which can be expressed as a sum of different 
multipoles. Nevertheless the order of angular scale of vari- 
ation for this function is similar to the dominant multipole 
probed. Only in close comparison with theoretical power 
spectrum estimates is an exact determination of the indi- 
vidual multipole terms of interest. For this we investigated 
the amount of "cross-talk" between the modes by varying 
(■max to the stability limits of the fit and studying the no- 
ticeable changes of the subtracted multipole modes, since 
the limitation to a range in £-space is the major problem 
for the fit technique. The changes to the results of our 
correlation analysis and dipole determination turn out to 
be insignificant at the present sensitivity level. In the case 
of fitting a dipole to a field the resulting parameters de- 
scribe the direction and amplitude of a local gradient and 
can only be compared for consistency with a whole sky 
dipole. In the case of the correlation analysis the errors 
induced by the effect are statistically taken into account 
by subjecting the simulated data to the same subtraction 
procedure. We find good agreement between the best fit 
multipoles and expected amplitudes from the correlation 
results on the region, or a subset of the region for which 
the fit is constrained. This, combined with the fact that 
the correlation function seems to be a fairly unbiased esti- 
mator for the amplitude of multipoles (Bunn et al. 1994), 
at least in the case of a power-law model for the CMB 
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fluctuations, gives us further confidence in the validity of 
the method used and we did not see the need to apply 
more sophisticated methods such as, e.g. constructing or- 
thogonalized functions on parts of the sky (Gorski 1994) 
or utilizing model assumptions in Wiener filtering tech- 
niques (Zaroubi et al. 1995), in the present state of the 
work. 



3.3. Quantifying Results 

For quantifying limits on the strength of the correlation 
(/?), we use a method which evaluates a likelihood distri- 
bution for (3, see e.g. Bennett et al. (1993). Assuming that 
the signal in the ROSAT data (X) and the COBE data 
(T) can be written in pixel space, leaving out the index i, 
as 

T = TcMB + T noise + P T \ X X 

(9) 

X = X X rb + X noise + Px\tT 

where Pt\x an d Px\t are the regression coefficients which 
couple the two maps by regarding one as template for the 
other. Clearly the COBE noise and CMB cosmic variance 
dominating the temperature distribution are the main 
source of confusion for the cross-correlation. The Poisson 
noise in ROSAT by photon number limitation is appar- 
ently a small confusion term. Shot noise is hard to dis- 
tinguish from the XRB since we are interested in a pos- 
sible correlation of ROSAT sources. From the results of 
testing with different source exclusion thresholds, we de- 
duce that shot noise is small after exclusion of a few very 
strong sources. A probably considerable confusion term 
for the cross-correlation is chance alignment of the struc- 
ture in the diffuse XRB, which we denoted by Xxrb (here 
meaning the uncorrelated part). Due to the complexity of 
the diffuse XRB, which consists of various typs of sources, 
more so than the CMB, no established model for the fluc- 
tuations exists. There is interesting work (e.g. Lahav et al. 
1996) to model the extragalactic XRB fluctuations, how- 
ever, these models have not yet been compared to data 
and hence are still somewhat uncertain. We believe that 
overall, the COBE noise and CMB cosmic variance cor- 
related with the real structure in the ROSAT maps are 
the dominant confusion terms, and for now, we have not 
attempted to model the XRB. Forming the correlation 
functions then yields the relations 

< XT > = < XTcmb > + < XT noise > + 

0t\ x <XX> 

(10) 

< TX > = < TXxrb > + < TX nolse > + 

Px\t <TT> 
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where < XT C mb >, < XT nolse >, < TX X rb > & 

< TX no i se > are assumed to be zero in the sense of a 
statistical average. This leads to the approximate relations 

Pt\x ~ < XT >o° / < XX >o° 

(11) 

0x\t ~ < XT > ° / < TT > o . 

For one particular realization this is only true within 
some error, which may be determined through Monte- 
Carlo simulations assuming models for Tcmb, T noise , 
Xxrb and X noise . We also take account of statistical bi- 
asing, which turned out to be small in comparison to the 
random errors. We assume now that Tcmb and T noise 
are the dominant sources of error compared to Poisson 
noise and random structure in the XRB, and approximate 

< T{X X RB + X notse ) > by < {Tcmb + T nolse )X >. Our 
method of determining (3t\x and the corresponding un- 
certainty is to minimize 

X 2 = ]T (< TX > -0r\ x < XX >)l 

kl 

(< TX > -f3 Tl x < XX >), (12) 

with M k i = < X(T C mb + ) >ki which can be de- 

termined from a distribution of realizations of the model. 
Assuming Gaussian errors, the probability distribution for 
/3 can be drawn from the % 2 -distribution as 

P(j3)dj3 cx e~^d0. (13) 
4. Results on Various Scales 

Since the aim of this analysis is to compare a CMB mea- 
surement to the diffuse XRB, some individual very strong 
X-ray point sources, which could influence even a statisti- 
cal comparison, are removed from the maps. Furthermore, 
the influence of different point source exclusion thresholds 
on the results were studied and found to be mostly in- 
significant. In a few cases, e.g. MK 421, a BL Lac object 
at z <~ 0.031, the removal of the source lead to a decrease 
of the correlation signal on the DMR beam scale. Com- 
parison with known strong radio point sources showed 
no significant contribution to DMR (Kogut et al. 1994), 
so chance alignment seems a likely cause. Nevertheless, 
systematic comparison between COBE data and candi- 
date radio-loud X-ray point sources seems a worthwhile 
check. The best fit DMR residual dipolc, which would in- 
troduce substantial correlation of no physical significance, 
had been removed, in addition to the standard removed 
dipole. 

4.1. Selected NGP Field 

Correlating the raw data (without excluding source 
contaminated pixels or subtracting structure) yields a 
marginally significant positive correlation on large angu- 
lar scales (figure 1), which appeared to be independent of 
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Subtraction of Gradients (DMR x ROSAT R5 & R6) 
70" < 1 < 250°, b > + 40" 
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20 40 
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Fig. 1. Cross-correlation function between the COBE DMR 
53+90 GHz and the ROSAT bands R6 (black) and R5 (white) 
in the selected NGP field (+40° < b, 70° < 1< 250°). The ef- 
fect of lowering the zero-lag amplitude by subtracting different 
multipoles related to a gradient on the field is shown (lowest 
panel). Subtracted are a best fit quadrupole from the Galac- 
tic custom cut DMR map, a best fit dipole which has Galac- 
tic signature from the Galactic cut ROSAT map (|b| > 20°) 
and a gradient as a dipole fitted on the field. The two latter 
both have Galactic signature (table 1). The 1-a error bands 
are taken from DMR noise + CMB simulations correlated to 
the ROSAT energy band R6. 



the following different procedures that had been applied. 
Different source exclusion thresholds in the ROSAT maps 
ranging from 0.3-1 cts s _1 were compared. After exclud- 
ing the strong source MKN 421 with 5.3 cts s _1 in the 
R6 energy band, the results were only marginally affected 
by different thresholds. Different sampling tests were un- 
dertaken, also showing stability of the result against 
small scale features such as point sources and noise. The 
maps were smoothed on various angular scales including 
smoothing of the ROSAT maps with the actual DMR 
beam (Kneissl & Smoot 1993), and Gaussian smoothing 
of both maps out to 20°, with the effect of smoothing 



the correlation function, but not significantly changing the 
correlated signal. 

The energy dependence in X-rays is found to increase 
from hard to soft energies. The frequency dependence in 
microwaves is somewhat unclear. There is a clear signal 
in both the 53 and 90 GHz channels and no signal at 31.5 
GHz. 

To determine the angular scale of the correlated signal, 
gradients were removed from the field. This was done in 
fitting dipoles onto the field in both maps and subtracting 
them. As a result the signal is reduced below the 1-a level. 
The multipoles on the sky dominating these gradients turn 
out to be of low order (figure 1). 

The gradient in the ROSAT selected NGP field has 
similar orientation as the whole map dipole (table 1), the 
positive pole lying near the Galactic center, and increas- 
ing in relative amplitude from hard to soft energies. From 
this energy dependence we derive spectral properties of 
the emitter that are in agreement with a 2 x 10 6 K equi- 
librium plasma, typical for Galactic emission. The energy 
dependence, however, is not compatible with the X-ray 
spectrum of the extragalactic XRB (Hasinger et al. 1993). 

The gradient in the COBE field is dominated by a 
quadrupole fitted to the COBE cut sky, a combination of 
the cosmic and the Galactic quadrupole. In the field, the 
cosmic quadrupole seems to dominate, which would not be 
inconsistent with the COBE frequency dependence of the 
correlated signal. A signal constant with frequency would 
be expected, but the 31.5 GHz channel could be confused 
by the increased noise level and Galactic contribution. 

Removing the gradients no significant correlation is 
left, and we can set upper limits of 4.5 /xK and 0.02 cts 
s _1 per pixel (95 % CL) in the R6 band on a correlation 
between the CMB and the extragalactic XRB on scales of 
7° - 40°. 



4.2. Cut Sky 

For completeness we present the correlation function on 
a large fraction of the sky, the DMR Galactic cut sky, al- 
though these results are not relevant for a comparison with 
the extragalactic XRB. We use the DMR cut sky for both 
data sets, since the most interesting question in this case 
is that of a correlation to Galactic features in the COBE 
DMR, which should be largely excluded by the custom 
cut (Banday et al. 1997), which is a straight |&| > 20° cut, 
with additional cutouts for Galactic structures found in 
correlation to the 140 ^m COBE DIRBE map ("flares of 
obscuration" in Scorpius, Ophiuchus, Taurus and Orion). 
We additionally experimented with |6| > 20° and 30°, and 
found, not surprisingly, a slight, not very significant ten- 
dency towards positive correlation with smaller cut angle. 
To minimize confusion, we subtract before the correlation 
in the cut sky case, in addition to the best fit DMR dipole, 
also the best fit ROSAT dipole. 



R. Kneissl et al.: Search for correlations between COBE DMR and ROSAT PSPC All-Sky survey data 



7 



Subtraction of the DMR Quadrupole (DMR 53 + 90 x ROSAT) 



DMR cut sky, R sources < 1 cts/s 




50 100 150 

Angle a ( °) 



Fig. 2. Cross-correlation between DMR 53+90 GHz and dif- 
ferent ROSAT energy bands including (black) and excluding 
(white) the DMR quadrupole. The DMR and ROSAT data are 
both taken from the same region, i.e. the DMR Galactic cut 
sky. Apart from the R6 band the narrow (solid line) and wide 
(dotted line) 1-a error bands are taken from DMR noise and 
CMB power spectrum simulations excluding and including the 
variance from the Q rm s-ps = 15.3 fiK respectively. In the R6 
band panel we demonstrate the 1-a error bands from (with in- 
creasing amplitude at zero-lag): DMR noise only, noise + CMB 
{£ > 2), noise + CMB (£ > 2) + Q rms (= 6 ^K) (dashed line) 
and noise + CMB (I > 1). There is no significant correlation, 
but a trend of increasing positive quadrupole correlation to- 
wards softer X-ray energy bands with about 68 % CL at R4 
against the cosmic variance of Qrms-ps = 15.3 fjK. The small 
excess in correlation on scales out to 10° most prominent in 
DMR x R6 has been found to be partly due to the LMC. 



The corresponding ROSAT sky is largely dominated 
by strong Galactic features such as, e.g. the Loop I with 
the prominent North Polar Spur, or the Eridanus enhance- 
ment, which are claimed not to be found in the DMR maps 
(Kogut et al. 1996) and therefore are sources of confusion. 
The limits we get confirm this view. The limits on the con- 
tribution of ROSAT to COBE are stronger compared to 
the ones derived on the selected NGP field, but this has to 



Subtraction of the DMR Quadrupole (DMR x R5) 
DMR cut, sky, R sources < 1 cts/s 



0.06 




50 100 150 

Angle a ( °) 



Fig. 3. Cross-correlation between ROSAT R5 and the individ- 
ual DMR frequency channels (GHz) including (black) and ex- 
cluding (white) the DMR quadrupoles. The DMR and ROSAT 
data are both taken from the same region, i.e. the DMR Galac- 
tic cut sky. The 1-a error bands are taken from DMR noise 
and CMB simulations including and excluding Q rm s~ps re- 
spectively. There is a trend of increasing quadrupole correla- 
tion towards lower DMR frequencies with a significant (> 95 
% CL) detection at 31.5 GHz. 

be taken with caution, since ROSAT is by a factor of <~ 2 
less constrained, due to increased uncorrelated (Galactic) 
structure in ROSAT. 

In figure 2 we find the DMR quadrupole to be posi- 
tively correlated to the ROSAT R4-6 templates. Here the 
outer 1-a error bands demonstrate the inclusion of the 
Qrms-ps — 15.3 /iK quadrupole. We stress here the fact 
that the huge cosmic variance of this term (~ 68%) intro- 
duces a large error, which is problematic, as the actually 
measured, uncorrected DMR Q rms is considerably smaller 
(15.6±5.4^K [31.5 GHz], 4.4±3.3mK [53 GHz], 0.0±3.0/uK 
[90 GHz]), because of an apparently anti-correlated align- 
ment of the CMB quadrupole (10.7 ± 3.6(random) ± 
7.1 (systematic) fiK) with the Galactic quadrupole (Kogut 
et al. 1996). For the energy dependence of the effect, we 
find agreement with a thermal 2 x 10 6 K spectrum and dis- 
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agreement with an extragalactic spectrum (a power-law 
with a photon index T w —2). 

Comparing the ROSAT R5 band to the different DMR 
frequencies (figure 3) we find a significant positive corre- 
lation (> 95 % CL) at 31.5 GHz due to a quadrupole term 
{QrS.F E = 10-3 ± 4 fiK, Q™? AT = 0.38 ±0.15 cts s" 1 per 
pixel). At the other frequencies a correlation due to the 
quadrupole is also positive, but not significant compared 
with the full CMB spectrum including the quadrupole. In- 
vestigating the spectral dependence we find the spectral 
index (i,T anten na oc v@ , within 1 a to be consistent with 

&ff ~ — 2.1 and p 'synchrotron ~ —2.8. 

When subtracting the best fit quadrupole, there is in 
no case any significant correlation left in the data. The 
limits are in the case of DMR 53+90 x ROSAT R6 on a 
30° cut sky 2.3 /jK and 0.037 cts s" 1 per pixel (95 % CL). 

From the systematics of the quadrupole correlation, we 
conclude that a correlation between Galactic quadrupoles 
in ROSAT and COBE seems likely, but which however, 
apart from the COBE DMR 31.5 GHz channels, can not 
be shown with sufficient significance against the cosmic 
variance of a quadrupole consistent with the CMB spec- 
trum. We see this result in agreement with the existence of 
a Galactic DMR quadrupole (Bennett et al. 1992; Kogut 
et al. 1996). 

4.3. Dipole 

The COBE dipole is the well known Doppler dipole due 
to the sun's motion with respect to the CMB rest frame. 
Any other residual dipole is an inseparable combination of 
imperfect Doppler dipole subtraction, Galactic dipole and 
CMB dipole. Thus a best fit dipole has to be removed from 
the maps, because the induced correlation has no physi- 
cal interpretation and thus confuses the results. Still the 
best fit dipoles in ROSAT can be determined and simply 
compared to amplitude and orientation of e.g. the DMR 
Doppler dipole, an expected Galactic dipole, etc. 

The dominant part of the correlated signal on the se- 
lected NGP field can be explained by a Galactic gradient 
in ROSAT. Comparing orientation and amplitude of this 
gradient expressed as a dipole on the sky, to dipoles fitted 
to other regions (table 1), we find it to lie in the general 
direction of the Galactic center and North Polar Spur re- 
gion and to have an amplitude of the dipole component 
towards the Galactic center (6n) which significantly in- 
creases towards softer X-rays, consistent with a Galactic 
energy spectrum. An increase in Galactic latitude towards 
harder X-rays in the selected NPG field, although also pre- 
ferred by fitting in the SGP field, has to be attributed 
to local Galactic phenomena, since it is not supported 
by a more global, joint fit to both fields combined. The 
X 2 /degrees of freedom (dof) indicates that the gradients / 
dipoles are clearly not the dominant structures, but they 
are fairly well defined in terms of the formal fit errors and 
insensitive to point source contributions (compare panels 3 



and 4 of table 1). For the determination of an extragalactic 
global dipole, systematic errors induced by Galactic fea- 
tures appear to be the overwhelming source of confusion. 
The dipole fits presented here can be compared to work 
by Freyberg et al. (1996), who investigated the ROSAT 
data with regards to a possible Galactic X-ray halo. 

5. Discussion 

5.1. Dipole 

Comparison of the dipole components is of interest, be- 
cause the firm detection of a cosmological X-ray dipole 
can either give confirmation to the velocity interpretation 
of the CMB dipole, or probe the distribution of matter 
on intermediate scales (100 - 1000 h _1 Mpc) and redshifts 
(z < 5). The velocity dipole amplitude can be calculated 
from the Compton-Getting effect (Compton & Getting 
1935) incorporating the relativistic effects of aberration 
and spectral shift to 

^ = (2 + r)^cos9. (14) 

The photon index (n(E) cx E~ r ) for the extragalactic 
XRB in the ROSAT energy bands is T ROS at = 2.0 ± 0.1 
(Hasinger et al. 1993) and the sun's velocity with re- 
spect to the CMB has been measured in the COBE 
DMR to be » = 369.0 ± 2.5 km/s in the direction 
(£,b) = (264°.31 ± 0°.17,+48°.05 ± 0.°10) (Lineweaver 
et al. 1996). With these numbers, the expected dipole 
rms amplitude, expressed as a percentage of the monopole 
term, is Dec = 0.08. This is more than an order of magni- 
tude lower than our typically measured dipole amplitudes 
(cf. to the last column of table 1, which gives the observed 
values expressed in the same manner). A dipole resulting 
from the distribution of matter is theoretically less well de- 
termined, but is assumed to have comparable amplitude 
to the Compton-Getting dipole (Lahav et al. 1996). Thus, 
the expected extragalactic dipoles have amplitudes lower 
than our observed values. 

Since we know from the energy dependence and orien- 
tation the Galactic origin of the fitted dipoles, our conclu- 
sion is that the ROSAT data over the energy range from 
0.5 - 2 keV are, in major regions of the sky, strongly dom- 
inated by a dipolar Galactic structure and a considerable 
effort would be needed to separate out an extragalactic 
dipole, if existent in the data. 

5.2. Galactic Quadrupole 

Strong Galactic correlation between the data should not 
be expected. A spatial correlation between regions of ra- 
dio free-free (Bremsstrahlung) emission and Hn regions 
(H Q ) is plausible. However, it is still unclear if the X-ray 
signature of Hn regions is dominated by emission from as- 
sociated hot gas or absorption by associated neutral gas. 
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R4 


4012 


36.8 


231.0±0.2 


-4.3±0.1 


5.9±0.1 


30.0±0.1 


351.8 


15.5 


8.7 


R5 


4012 


54.3 


275.7±0.2 


-5.3±0.1 


4.2±0.1 


26.4±0.1 


348.7 


12.4 


7.7 


R6 


4012 


34.8 


315.9±0.2 


-4.2±0.1 


2.5±0.1 


16.9±0.1 


346.1 


11.6 


5.0 


R7 


4012 


7.1 


169.3±0.2 


-3.5±0.1 


0.6±0.1 


6.0±0.1 


329.4 


7.2 


2.0 



Table 1. Dipole fits to various subsets of the ROSAT PSPC All-Sky Survey data. The fields NGP (b > +40°) and SGP (b 
< —40°) are in the longitude ranges indicated. The source exclusion thresholds (sre) are in cts s _1 . h im are the coefficients 
of the real valued spherical harmonics, whereas b\ m is 100 x b\ m / boo, the dipole coefficients in percentage of the monopole 
term. (I 11 , b 11 ) gives the best fit position for the positive pole in Galactic coordinates and D 2 = 1/(471") y^ m __j b\ m is the 
amplitude. Although we show the (small) fit errors on the 6i m resulting from ROSAT noise only, we do not translate them into 
the coordinate values as we see the results dominated by systematic errors. For details see text. 



Only in the first case would a positive correlation be ex- 
pected. 

Supernova remnants are prominent Galactic features 
in the soft X-rays, and enrich the Galactic medium with 
compressed magnetic fields that are responsible for syn- 
chrotron emission also at microwave frequencies. However, 
under the usual assumptions regarding the spectral be- 
haviour of these synchrotron sources, a strong contribu- 
tion at COBE frequencies is not to be expected. This is in 
agreement with the lack of correlation (Kogut et al. 1996) 
between COBE DMR and the template for synchrotron 
radiation, the 408 MHz map (Haslam et al. 1982). The 
North Polar Spur, which traces the rim of the nearby 



superbubble Loop I (e.g. Egger et al. 1996), is however 
the most prominent synchrotron source and one positive 
pole of the correlated quadrupole lies in the correspond- 
ing direction. The amplitude of the correlated signal in the 
COBE DMR is not incompatible with the signal expected 
from this object in the microwave band, given the uncer- 
tainty in the spectral index. Since the North Polar Spur 
is also known as a very bright X-ray source it is possible 
that this feature is associated with the correlation. 

Although the contribu- 

tion of the correlated quadrupole to the CMB channels 
in the DMR is comparatively small, it will be interesting 
to investigate the meaning of the detection at 31.5 GHz 
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EFFECT / 
SOURCE 


SIGN 


ANGULAR 
SCALE 


FREQU. DEP. 
microwave 


ENERGY DEP. 
X-ray 


AUTHOR 


Galaxy 

geometrical ? 


+ 

(-) X-abs. 


large 


ftsynch, f f ,dust ^ 

{-2.8,-2.1,1.5} 


thermal 
0.3 keV 




SZ thermal 
(clusters / super- ) 


(Rayleigh-Jeans) 


(< 10' clust./ 
< 5° c-corr.) 


y - distortion 


thermal 
10 keV 




SZ thermal 
local group halo 


(Rayleigh-Jeans) 


large 


y - distortion 


thermal 
1 keV 


[1] 


X-ray/radio 
point sources 


+ 


small 


fiat, a < 0.5 
(10 - 100 GHz) 


T ?a 2.0 ±0.1 

I = Io(E/E Q )- r 


[2,3,4] 


ISW / RS 

in A / open universe 


+ 


large 
I « 10 


Planck 


r « 2.0 ±o.i 

/ = I Q (E/E )- r 


[5,6] 



Table 2. Overview of different effects introducing possible correlations between microwave and X-ray data. Authors: 
[1] Suto et al. 1996, [2] Franceschini et al. 1989, [3] Franceschini 1995, [4] Laurent-Muehleisen et al. 1996, [5] Crittenden 
& Turok 1996, [6] Kamionkowski 1996. 



further, by comparison with the Galactic synchrotron ra- 
diation template and the template for Bremsstrahlung and 
dust radiation (which is the COBE DIRBE 140 /im map). 
These topics will be addressed in a forthcoming paper. 
By doing this, a separation between synchrotron radia- 
tion and Bremsstrahlung might be possible. The answer 
to the question of whether a physical effect or purely geo- 
metrical alignment is responsible for the correlations could 
be decided. 



5. 3. Derivation of Limits 

On the selected NGP field some indication of a positive 
correlation of large angular scale was found. Since the cor- 
related gradient in ROSAT is of Galactic origin, chance 
alignment with the COBE CMB quadrupole, and pos- 
sibly some spatial alignment with the COBE Galactic 
quadrupole (see section 5.2) turns out to be the cause, 
and not an extragalactic correlation. After subtracting 
the dipole and quadrupole terms from a Galactic cut sky, 
or related large angular scale gradients from the selected 
NGP field, no significant correlations are found in the 
data. This will enables us to set limits on a number of 
possible mechanisms for introducing correlations as de- 
scribed in the introduction and summarized in table 2. We 
note here, that all our findings are in no conflict with the 



standard cosmological interpretation of the COBE DMR 
measurements. 

The large-scale, correlated feature we are seeing, the 
presumable Galactic quadrupole, cannot originate from 
the SZ effect, as would be the case for the Local Group 
halo, since the observed correlation is positive. So at the 
quoted sensitivity, no Local Group halo is found. By mak- 
ing use of the specific halo template the sensitivity might 
be improved. 

From the temperature limits on the selected NGP area 
we can infer limits on the Comptonization parameter y on 
scales of 7°-40° from distortions in the Rayleigh-Jeans 
regime 

6T/T = -26y (15) 

and find 5y < 8 x 10~ 7 . More specific constrains will be 
derived by combining the DMR and ROSAT limits with 
model assumptions about the distribution of the gas and 
its properties. For the gas halo model adopted by Soltan 
et al. (1996a) as one explanation for the extended corre- 
lated X-ray emission around clusters of galaxies, we can 
limit the assumed temperature of smoothly distributed 
gas to less than 2 keV (95 % CL). Miyaji et al. (1996), 
when correlating ROSAT to HEAO, found excess fluctu- 
ations at ROSAT energies in comparison to predictions 
from the population synthesis model by Comastri et al. 
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(1995). These fluctuations can be explained with excess 
emission by the gas halo model only when a temperature 
of at least 2 keV is assumed. Combined with our limits, 
this leaves a very narrow, possible temperature range for 
the gas. 

Also the combined contribution to the microwave and 
X-ray fluxes of a population of radio-loud AGNs can be 
limited, taking into account that the sensitivity for the 
correlated locations is limited to the zero-lag convolved 
with the DMR beam. 

No indications for a correlation through the ISW have 
been seen, which should dominate in the high DMR fre- 
quencies and hard X-rays. In this case, where we arc in- 
terested in correlations between the CMB and the ex- 
tragalactic XRB at a range of scales from a few degrees 
up to the quadrupolc, another search strategy could be 
adopted, correlating the selected NGP field to the DMR 
cut sky with the effect of slightly improving the sensitivity 
on scales of the size of the field, as tests have shown, and 
quantifying the sensitivity on scales above. Since we do not 
know how well the fluctuations in the ROSAT XRB trace 
the projected gravitational potential, although rough as- 
sumptions can be made from the properties of the resolved 
sources (Comastri et al. 1995), we do not give any limits 
on a A-universe in this analysis. 

Comparing the limits derived on the selected NGP 
field and on the 30° cut sky, we find them to differ 
only marginally. Two effects, reduction of v^-noise in 
COBE and increase of predominantly uncorrelated Galac- 
tic structure in ROSAT, when going from the smaller to 
the larger area, can be seen. For comparison with the ex- 
tragalactic, diffuse XRB the limits from the selected NGP 
field are preferred. 

We qualitatively discussed how the limits we found 
constrain possible correlation mechanisms. A more quan- 
titative study should make use of specific predictions by 
each individual effect for e.g. frequency and energy depen- 
dence, angular scale, spatial orientation etc. which could 
result in stronger constraints. This analysis is in progress 
and will be presented in a future publication. 

On the observational side, great improvements for this 
kind of analysis can be expected from future satellite 
experiments such as the CMB missions MAP and CO- 
BRAS/SAMBA (Bersanelli et al. 1996), compared to high 
resolution X-ray observations such as XMM and AXAF, 
and the hard X-ray survey ABRIXAS (Fricdrich et al. 
1996). In terms of statistical analysis a comparison be- 
tween COBRAS/SAMBA and ABRIXAS appears most 
promising, since ABRIXAS will be the most sensitive all- 
sky X-ray survey (0.3-10 keV) and COBRAS/SAMBA 
will be the CMB mission with highest angular resolu- 
tion and best frequency coverage, enabling a separation 
between different frequency dependences such as Planck 
spectrum, y-distortions and radio source spectra to be 
achieved. 
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